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NOMENCLATURE 


average wall superheat (K) 

temperature at the base of the fins (°C) 
fluid saturation temperature (°C) 
average wall temperature (°C) 


heat-transfer rate from the cartridge heater 
(W) 


heat-transfer rate through unenhanced end (W) 
total heat-loss rate (W) 

heat flux from the Gewa-T surface (w/in 
average heat-transfer coefficient C/ ae 
thermal conductivity of the copper (W/m.K) 


thermal conductivity of Freon-ll3)ei¢u 
(W/m. K) 


volumetric thermal expansion coefficient 


(1/K) 

kinematic viscosity ae 
thermal diffusivity ee 
density of Freon-113 (kg/m?) 
Nussel number 

Rayleigh number 

Prandtl number 


voltage across the cartridge heater element 
(volts) 


voltage drop across the precision resistor 
(volts) 


current through the heater element (amps) 


me 


resistance of the precision resistor (ohms) 
outside diameter of the Gewa-T tube (m) 


diameter at the position of the thermocouples 
(m) 


€Giametervat the base of the fins (m) 
outside diameter of the unenhanced ends (m) 
inside diameter of the unenhanced ends (m) 


outside wall perimeter of the unenhanced 
ends (m) 


fin height (m) 
length of enhanced test section (m) 
length of unenhanced ends (m) 
corrected length of unenhanced ends (m) 
tube thickness of unenhanced ends (m) 

2 


surface area at the base of the fins (m ) 


cross-sectional arga Qu cube wall atk the 
unenhanced ends (m*) 


outside surface area of the unenhanced ends 
(m<) 


gravitational acceleration (m/s) 


tial 


ACKNOWLEDGMENTS 


Above all, the author must thank for the continuous 
encouragement he received from his wife. 

For their infinite patience and accurate guidance 
throughout this investigation, the author is sincerely 
grateful to Dr. Paul Marto and Dr. A.S. Wanniarachchi. 

For their cooperation in furnishing samples of Gewaog 
tubing, the author wishes to thank Wieland Werke AG, West 


Germany. 


eZ 


ieee RODUGCT LON 


A. BACKGROUND 

Lieweneboy eCbisis tena began im 1974 has made a consider- 
able impact on the design and operation of energy-conversion 
devices. Today, designers are challenged to design more 
and more efficient heat-transfer systems. Among the numer- 
ous areas that designers are faced with this challenge, 
boilers or evaporators used in many engineering applications 
demand considerable attention. The use of advanced boiling 
surfaces may lead to significant improvements in efficiency 
and/for significant savings in the initial cost and equipment 
sizes. Also, due to the ever-decreasing sizes of electronic 
Semponencs, the presence of very high heat fluxes has made 
the cooling of such components a challenging problem. Marto 
and Lepere [Ref. 1] stated that the use of advanced boiling 
Surfaces may have potential promise in the field of electronic 
cooling. However, the lack of generalized performance data 
and useful theoretical models for such advanced surfaces 
have still kept designers away from the utmost use of these 
Surfaces. 

Due to a lack of complete understanding, the theoretical 
treatment of the nucleate pool-boiling process, even ona 
smooth surface, is extremely complex. The existence of a 


large number of semi-empirical formulas is probably a very 


io 


good indication of the complex nature of the nucleate boil- 
ing process. Chongrungreon and Sauer [Ref. 2] summarized 
nine(9) correlations for nucleate pool boiling 2 5si7g. 
component fluids from smooth tubes. They stated the inade- 
quacy of any of these correlations when compared with 
experimental data. 

Nucleate pool boiling is strongly dependent on the Sums 
face structure; i.e., the heat flux increases with increase 
in stable nucleation site density. The size, shape and den- 
sity of such nucleation sites are almost impossible to pre- 
dict. However, the entire nucleate pool-boiling process is 
dependent on the performance at these nucleation sites. It 
is well known that the ability of these sites to develop and 
trap small vapor bubbles is the reason for improved performance 
in nucleate pool boiling in comparison with natural-convective 
De anc. | 

When a bubble is being generated at a surface nucleation 
Site (or a pore), the wall superheat is given by the follow- 


ing approximate theoretical equation [Ref. 3]: 


ZG To 
Oye 
Thus, for a given fluid, the wall superneaeyecanes eae 
creased (boiler performance increases with decrease in wall 
sSuperheat) with increase in pore radius. Of course, this 
radius must have a practical upper limit since too large a 


pore would not trap a vapor bubble. 
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Nishikawa and Ito [Ref. 4] recently discussed two methods 
memaugment nucleate boiling: (1) use teflon-coated pits to 
reduce wettability of boiling ligud, and (2) manufacture 
surfaces with numerous re-entrant cavities that have the 
emi tO trap vapor and keep the nucleation sites active. 
Due to the high wettability of low surface-tension fluids 
(such as fluorcarbons), the use of teflon-coated pits has not 
been successful. Thus, the surfaces with re-entrant-type 
cavities are the viable candidates for low surface-tension 
rcquids . 

Taking advantage of the concept of re-entrant cavities, 

a number of enhanced boiling surfaces have been invented 
Guring the last two decades. They include the Linde High 
Flux surface of Union Carbide, Gewa-T surface of Wieland 

and Thermoexcel-E surface of Hitachi. All of these surfaces 
are provided with re-entrant cavities through specialized 
manufacturing processes. These re-entrant cavities introduce 
Temmmore parameters to be desired. 

In view of the inadequate understanding of the nucleate 
pool boiling process on smooth tubes, a thorough understand- 
ing of the boiling process on advanced surfaces may take 
many more years to come. Thus, present-day designers have 
metry little information when considering advanced boiling 
Surfaces. It 1S very important that researchers gather a 
vast amount of experimental data covering a wide range of 


heat-transfer fluids. Such extensive data will enable the 


ies 


development of generalized correlations or successful 
theoretical models. 

Yilmaz, et al. [Ref. 5] compared the nucleate pool-bomiame 
heat-transfer performance of three copper tubes: Lindewiaa 
Flux, Wieland Gewa-T and Hitachi Thermoexcel-E. Their data 
taken with p-xylene showed that while all tubes performed 
better than smooth tubes, the best performance was with the 
Linde High Flux surface. The Thermoexcel-E tube performed 
better than the Gewa-T tube at low heat flux and Dboeh oem 
formed similarly at high heat flux. These data were taken 
under the fully-established boiling regime; thus, no intone 
mation was provided on the hysteresis behavior or temperature 
overshoot of these surfaces. 

Bergles and Chyu [Ref. 6] compared the nucleate pool- 
boiling heat-transfer performance of four different Linde 
High Flux tubes to a plain copper tube in distilled water saw. 
in R-ll3. They used three different methods of aging the test 
surface prior to collecting data. They reported significant 
temperature overhsoots of both plain and High Flux surfaces 
in R-1l13. They also reported that these overshoots were 
sensitive to aging, initial subcooling and rate ereee 
imerecase. 

Marto and Lepere [Ref. 1] obtained results on the heat- 
transfer performance of three heat-transfer surfaces: Linde 
High Flux, Thermoexcel-E, Gewa-T, and a plain copper sur- 


face in R-1l13 and FC-72. They reported a two to tenfola 
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increase in the heat-transfer coefficient when compared to 
plain surface. The High Flux surface was most effective over 
PmeGedderange Of Neat Eluxes. The Thermoexcel-E surface 
showed similar gains in heat-transfer coefficient to that 
Se the High Flux surface below 10 kW/m° . The Gewa-T surface 
was not as effective as the othér surfaces at low heat fluxes 
pms performed Comparably at high heat fluxes (near 100 kW/m“) . 
Marto and Hernandez [Ref. 7] studied the boiling per- 
formance of a Gewa-T finned surface in Freon-113. They 
examined the influence of placing aluminum shrouds around the 
test tubes, and found that the thermal performance of these 
surfaces were dependent upon the liquid-vapor flow within 
the channels between neighboring T-shaped fins. They also 
examined the influence of the fins by progressively machining 
away the fin height to arrive at a smooth cylindrical surface. 
They reported that the addition of the "T-caps" to straight 
fins produced the most significant improvement in heat- 
transfer performance when compared to a smooth tube. The 
use of shrouds increased the heat-transfer coefficient at 
all practical heat fluxes, but decreased the heat-transfer 
coefficient at very high heat fluxes when compared to the 


unshrouded surface. 


Pee ©HESIS OBJECTIVES 
Easeowoneeie foregoing discussion, the major motivation 
in this work was to further study the boiling heat-transfer 


performance of Gewa-T surfaces in R-113, and to attempt 


ee 


to improve this performance through the use of further 
enhancing techniques. 
There were four major objectives in this thesis study: 

1. Take data to verify results of Hernandez [Ref. 8] on 
a standard 21.2 mm 0.D. Gewa-T tube with and without shrouds. 

2. Use wire wraps in the inter-fin area to improve 
nucleate-boiling heat-transfer performance from Gewa-T 
tubes. 

3. Take data using the combination of best shroud (foam 
in 1) and the best number of wire wraps (found in 2) on 
Gewa-T tubes. 

4. Take data to study the effect of fin-tip spacing and 
Fin density on the heat-transfer performance of Gewa-T 


tubes. 


ike 


II. EXPERIMENTAL APPARATUS 


Pe REQUIREMENTS GOVERNING DESIGN 
The purpose of this research was the observation of 

nucleate pool boiling from different Gewa-T surfaces. In 
Order to get a more comprehensive knowledge about the physi- 
cal mechanisms that govern this phenomenon, it is necessary 
to measure the parameters involved which affect the nucleate 
boiling process. The selected parameters to be obtained 
are the following: 

1. Barometric pressure 

2. Test section wall temperature 

eee Flugad bulk temperature 

we vapor temperature 


oa Surface heat flux 


fee LEST APPARATUS 

The experimental apparatus designed by Lepere [Ref. 9], 
and modified by Hernandez [Ref. 8], with further minor modi- 
Precations, was used during this thesis work. To gather data 
more conveniently at a faster rate, a Hewlett Packard model 
9826 computer and an HP-3054A automatic data-acquisition 
system were added to the experimental facility. Lepere 
[Ref. 9] and Hernandez [Ref. 8] provide detailed descriptions 
Of the apparatus; therefore, only a brief description is 


presented in this thesis. 


ily, 


Figure 1 shows a schematic of the test ™apparatus ee 
boiling surfaces consist of a cylindrical Copper usen. een 
a smooth interior surface and a commercially-available T- 
shaped finned surface manufactured by Wieland-Werke, AG. Each 
tube was conditioned for the insertion of thermocouples by 
drilling eight longitudinal holes as indicated in Figures. 
and introducing the respective copper-constantan thermocouples 
centered by means of capillary copper tubes. In the center 
of the tube a cartridge heater was placed to act as the heat 
source as shown in Figure 3. 

Both ends of the test section were machined to remove 
the enhanced surface on the outside and to increase the 
interior diameter. This reduced the cross sectional area 
and thereby minimized the heat transfer from the ends of the 
test section. A Teflon plug was inserted into the free end 
for sealing and heat transfer reduction purposes. The other 
end of the test section was connected to a long copper ex- 
tension tube through a 90-degree elbow, which allowed the test 
section to be placed in a horizontal position within a thick- 
walled Pyrex glass vessel. Since one of the requirements 
of the apparatus was to maintain the liquid at saturation 
temperature and atmospheric pressure, a plate heater and a 
vent tube were provided to satisfy these conditions. In 
addition avoiding excessive loss of vapor, a water-cooled 
copper coil was placed in the vapor space of the vessel, to 


act as a primary condenser. The vessel was covered with a 
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Schematic of Test Apparatus 
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Plexiglas lid, which had been fitted with a tubberso--aaee 
Several holes were drilled in the lid to fix in place the 
test section extension, vent tube, connections fOr prime 
and secondary condensers, and two more thermocouples to sense 
both the bulk liquid and vapor temperatures. The Pyrex 

glass vessel with the attached fittings was placed on the 
plate heater, which served to maintain the liquid at satura- 
tion temperature. A secondary condenser was mounted external 


to the vessel. 


Cr INSTRUMENTATION 

The power supply was connected to the cartridge heater 
through a voltage regulator and then to a Variac scaled from 
O to 120 VAC. This arrangement allowed fine adjustment for 
a stable feed voltage during the runs. A precision resistor 
of 2.031 ohms was connected in series with the cartridge 
heater to provide a méasure of the Current Elow by Teaqmag 
the voltage drop across this element. A digital voltmeter 
was used to measure both the voltage across the resistor and 
the voltage across the cartridge heater. All the thermo- 
couples were connected to a computer data acquisition system, 
and ambient conditions and voltage readings were supplied 
manually. The microcomputer processed all this information 


through a program included in Appendix A. 


D. CALIBRATION OF THERMCGSUELES 
A total of fifteen wall thermocouples in tubes 1 and? 


(one thermocouple in tube 1 was defective, hence it was 
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excluded), and the two thermocouples for vapor and liquid 
temperature measurements were calibrated in an insulated, 
Stainless-steel bath. The tubes were completely immersed 
into the bath, which was provided with a motor-driven mixer. 
The bath temperature was measured using a platinum-resistance 
thermometer accurate to +0.01 K. 

Since the temperature measurements during this investiga- 
tion ranged from 45°C to 60°C, thermocouple calibration was 
performed in the range of 45°C to 70°C. To begin the cali- 
bration process, hot water (at about 75°C) was added to the 
masulated container. Following proper mixing, all thermocouple 
readings were automatically recorded through the data- 
acquisition system. The bath temperature was measured using 
the platinum-resistance thermometer. To obtain other lower 
temperatures, small quantities of water were added to the 
bath, allowing 5 minutes for temperature stabilization of the 
test section. 

The discrepancy (1i.e., the platinum-resistance thermometer 
reading minus the thermocouple reading) was plotted against 
the thermocouple reading as shown in Figure 4. The second 
order polynomial curve shown in this figure was generated 
using the data from all seventeen thermocouples. However, 
for clarity, only three thermocouples were selected for data 
points from each tube. It is evident that the thermocouples 
in both tubes show similar discrepancies, the thermocouple 


readings have only about +0.1 K discrepancy around the 
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calibration line. Based on this fact, no calibration was 
performed for the other five tubes. Instead, the calibration 


curve generated for tubes 1 and 7 was assumed for all tubes. 
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Iii. EXPERIMENTAL PROCERGRIS 


A. NORMAL OPERATION 

The Procedure B employed by Marto and Hernandez [Ref. 7] 
has been demonstrated to be the most adequate to carry out 
the boiling runs. Because of the large influence of surface 
past history on boiling incipience, this procedure was estab- 
lished to allow repeatability and permit comparison of the 
obtained data for different enhancements. 

The procedure is as follows: 

1. Once the test section has been immersed in the pool, 
it was subjected to one hour of pre-boiling by setcingm@iane 
cartridge heater to give a heat flux of about 30 kW/m* . The 
plate heater was adjusted to the maximum temperature for this 
one-hour period. This initial vigorous boiling served to 
degas the liquid and to force the noncondensable gases out to 
the atmosphere. 

2. After the initial aging process, both test-section heater 
and plate heater were secured, allowing cooling of the liquid 
and the test section for about 30 minutes. 

3. The plate-heater voltage was then adjusted in order to 
maintain the liquid at saturated conditions and the cartridge 
heater was turned on to the first setting of 12 volts. 

4. For all the consecutive settings during both increasing 
and decreasing heat fluxes, the boiling was allowed to 


Stabilize for five minutes at each power setting. 
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5. At a given power setting, the following data were 
recorded: heater voltage, precision-resistor voltage, vapor 
BenperatlmemilagquIcd outk temperature and wall temperatures 


me the teSe section. 


jee GEWA-T SURFACES 

Three groups of surfaces were tested. The first. surface 
was the same as used by Hernandez [Ref. 8] with a fin density 
of 0.75 fins/mm, an outside diameter of 21.2 mm, and an 
inter-~fin gap of 0.25 mm. The second group of three tubes 
had 0.75 fins/mm, an outside diameter of 23.1 mm and inter- 
Ppuecgaps Of 0.15, 0.25 and 0.35, and the third group of three 
tubes had 1.02 fins/mm and all other specifications were the 
eame aS the second group: 

The geometry of the enhancement of these seven Gewa-T 
surfaces was observed by making a longitudinal cut of the 
tubes and photographing them with 50X magnification (Figure 5). 
As can be seen, the tubes with 0.75 fins/mm had an approxi- 
mately circular cavity shape, while the 1.02 fins/mm tubes 
had a somewhat elongated cavity shape. This difference in 
Shape was probably due to the technique used in construction 
to obtain the various specifications. All the T-shaped fins 
mea a fin height of approximately 0.95 mm. 

The exact determination of diameters at the base of the 
fins was carried out from the magnified photographs of the 
cavities. This diameter was determined by subtracting twice 


the fin height from the measured outside diameter: 
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(b)  Gewa-T Tube Number 5 


Figure 5. Optical Micrograph of Gewa-T Fin Profile, 50X 
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mesulting in the dimensions presented in Table l. 


Peek 


Dimensions of Gewa-T Finned Tubes 


TOBE FINS/ PL Nese PNT tee GH 1 ee Ds D> 
INO’. mm GAP (mm) (mm) (mm) (mm ) 
i ORS G25 02865 Zee 19 347 
i Oats Oe Oo c0 24.99 Zoe 5 
8 O75 Ot Or 60 2a oo Do O 
4 ORS Oe) Oe 18 Care 237305 
> iO ees i050 2750 2a 20 
6 02 O25 10 50 meee 237.20 
i OZ roo ie 25 Dao) 2305 


For simplicity, the diameter at the base (D.) was assumed 


memoe 19.5 mm for Tube 1 and 23.1mm for all other tubes. 


fe oEAT—-FLUX CALCULATION 

meecording to earlier description about test-section con- 
struction, the cartridge heater acts as the heat source, 
inserted into the tube, and variations in heat-transfer rate 


Becemade througn different power settings with a variac. 


3k 


where: 


Q = heat-transfer rate from the cartridge 
H 
heater (W) 
V = voltage across the cartridge heater 


element (volts) 


i = current through the heater element (amps) 


This circuit is connected in series with a precision 
resistor, which provides the means to compute the electrical 


current by the measurement of the voltage drop across the 


resistor: 
Vee ® 2) 
where: 
V = voltage drop across the precision 
Res 
resistor (volts) 
R = resistance of the precision resistor (ohms) 


= 2,03). 20hnms 


Both ends of the Gewa-T tube were machined to provide a 
smooth surface and a reduced wall thickness in order to mini- 
mize longitudinal heat conduction and therefore the heat 
losses from the unenhanced ends. The geometry of these ends 


is shown in Figure 6, where: 


uy 
il 


base of the T-shaped fins 


T = temperature at the base of T-shaped fins 
LC} 


B)/7 


D 
3 ee 


a Se 


Figure 6. Geometry of Unenhanced Ends 


Meeechis cylindrical geometry: 


Si ee 2 
Aa = 4 (dD. D.) (3) 
A = PXL (4) 
where: 
AG = cross-sectional area of tube wall Cm) 
Da = tube outside diameter (m) 
D. = tube inside diameter (m) 
Ag = EmoowOltoidemoUEnace area a 
P = tube outside wall perimeter (m) 


It was assumed that the temperature at the base of the 
thin-walled fin was equal to the temperature at the base of 
the T-shaped fins. Therefore: 
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O(x) = rt) en (5) 


SAT 
In general, 
Z 
os - Gre = 0 (6) 
@ pe 
where: 
2 h P 
m = See (7) 
KAq 
Therefore 
_ 2 fees ee ces 
Qn = QO = ~KAG ax — TRAC ax ‘4 
x=0 x=0 


Now, assuming negligible convective heat loss from the 


thane We aeaso: 


de 


x= 
and 
6 . “€osh mi ba~ 
ag = Cosh mL (10) 


For this case 


O = Vh Pk A So manns mL, ( 115) 
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where: 


J| 


Assuming the 
of the eight cal 


measurements, le 


(12) 


heat-transfer rate through unenhanced 
end (W) 


fluid saturation temperature (°C) 


thermal conductivity of the copper (W/m.k) 


average heat-transfer coefficient (W/m? . K) 


wall temperature is the same as the average 
ibrated wall temperature thermocouple 


ss temperature drop due to radial conduction 


from thermocouple position to the base of the T-shaped fin, 


In(D,/D,) 
Sow aver oH 270K, ae! 
8 
T —e), Da) 7820 (14) 
AVG aq; 
where: 
Dy = diameter at the position of thermocouples 
(m) 
Tayc = average wall temperature (°C) 
For the geometry shown in Figure 6: 


(ales?) 


& =D. = oF (Gp 


O a 
where: 
La = corrected length (m) 
Li = length of unenhanced ends (m) 
it = tube thickness (m) 


The average difference between the wall temperature and 
the saturation temperature may now be determined from the 


following equation: 


a 1 a Cosh m(L,-x) 

6 = =— — de (176 
La 0 Cosh mL 

6 = —— Tanh mL ( ee 
L C 

ee 
9 = Ty - Tsar so 
where: 
(To-Toap) = average difference between surface 


wall temperature and fluid 
Saturation temperature (kK) 


Since the unenhanced ends are very long and the thickness 
is very small, only very little heat transfer would take 


place at these ends. Thus, these ends undergo free convection, 
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Pore vaten tne Churchill and Chu correlation as stated in 


[Ref. 10] for the average Nusselt number was assumed: 





m7 6 
O2eas7 pbs 5 
Nu = 0. 60 +... —-— == 
Dy fl ecomaeey ery 294 °7 7 | 
ome Pee ee 
Oo 
The average Nusselt number is: 
_ h Do 
Nu, = e 
Oo 
where: 
Kk = thermal conductivity of Freon-113 
pelving for h: 
0-388. Ral/® 
= ¥ 60+ “0 |? 
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(W/m. K) 


As) 


where: 


Pr = Prandelenuneer 
Vv = kinematic viscosity (ayes 
1 = thermal diffusivity (m*/s) 
Ray = average hay lergn nung 
O 
g = gravitational acceleration (m/s*) 
B = volumetric thermal expansion coefficient 
(ye) 
9 = density of freon-113 (kg/m?) 


Now, substitution of Equations 7, 18 and 24 into EQuatio@nmee 


results. ins 


— 1/2 IAG 
3 RiP a 
g8D°9,Tanh{ KAQ ‘al 2 
Ones7 - RP ) 1/2 (2G 
C'KAG J 
— k 
h = = <0.60 + 
D6 il + (0.559/7en) = eee 
Equation 26 is solved for h by iterative tecnmuame 

within a range of precision of 0.001. Knowing the value of 


natural-convective heat-transfer coefficient at unenhanced 
ends, the total heat-loss rate is calculated from Equation 
ll, and the heat-transfer rate through the enhanced Gewa-T 
surface is obtained by subtracting the total heat—lossiiage 


from the total heat-transfer rate: 


Qross 205 (27) 
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Op MO 205 (28) 


Hinatty, ene heat flux from the Gewa-T surface at the 





base of the T-shaped fins is: 


g _ Q/Az (29) 


A, eae (3.03) 


Be, 


IV. RESULTS AND DISCUSS Wer 


A. DESIGNATION OF DATA FILES 

A total of seven Gewa-T tubes were tested during this 
investigation. As can be seen from Table 1, the primary 
geometric variables were: fin-tip gap, fin @demoue aa. 
outside diameter of tubes including fins. 


The general designation for data files is: 


GTNPnn 


where: 


GT indicates Gewa-T surface 


N specifies the number of the tube tested 


(Table 1) 

P indicates the position of the test section 
(1.e., H for horizontal and V for vertieane 
and 


mn indicates the corresponding cumulative run 
number. 

For example, GT5H28 specifies tube number five tested in the 
horizontal position during run number 28. From Tastee 
the corresponding specifications for this test section are: 
25.3 mm outside diameter, 23.1 mm diameter at the base of the 
fins, 1.02 fins/mm and inter-fin gap of 0.15 mm. While the 
file name just presented contain data in the very raw form, 


the processed data (i.e., q versus AT) were stored ina 
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Separate file. This plot data file could easily be recognized 
bye tne “P" added at the beginning of the raw data file name. 
In this manner, the plot data file corresponding to the above 
raw data file has the name PGT5H28. 

PeemOQtalwoOn thirty Gqitrerent data runs were Carried out 
during this investigation. These runs break down in the 


following manner: 


Run Number Tube Number lestecondi tion 
1 - 3 al Plain —- Vertical position 

4 - 10 Hl Plain - Horizontal 
ees 2 1 dk With shrouds and/or Wires 

Zea 2 3 2 Pee ei ZOn tc ll 

24 = 25 3 Plot etor Zonta ll 

26 4 Pilati iera Zonta | 

Za 2 2 Plaeoe—e HOD Zonta L 

29 6 PILeivol — |slebentyeseusr- wk 

S10 7 Poti ior LZOnta L 


eee CLRCUMPERENTIAL UNIFORMITY OF HEAT FLUX AROUND HEATER 
ELEMENT 


During the early stages of this investigation, it was 
thought necessary to test if the heater element would provide 
a non-uniform heat flux along the circumferential direction. 
Such a non-uniform heat-flux variation was assumed possible 
based on the heater-element construction. Therefore, tube l 


was immersed vertically in the R-113 liquid pool and three 


41 


runs were made. All of these runs showed only a negligible 
variation of the wall thermocouple temperatures; maximum 
variation was +0.25 K. These runs clearly showed that the 
heater element indeed provided a uniform heat flux alongmaa. 
circumferential direction. No detailed data are presented 
in this thesis for the tube in the vertical position. All 
other results presented herein are for the tubes in the 
hori Zonta Psoestetene 
C. LONGITUDINAL UNIFORMITY OF HEAT FLUX ALONG THE TES? 

SEGLEION 

During most of the data runs with the tubes in the hori- 
zontal position, somewhat unexplainable phenomena were 
observed. As the heat flux was gradually increased, bubble 
formation began simultaneously at the two ends, while the 
center portion of the tube underwent only natural convection. 
Approximately 40 percent of maximum heat flux was necessary 
to obtain nucleate boiling from the entire tube. Since all 
the wall thermocouples are located at midway along the tube, 
the average wall temperature values provided in all data 


runs may be in slight error. 


D. PLOT ANALYSIS FOR PLAIN SURFACE 

Figure 7 represents a typical curve for nucleate peat. 
boiling process of Gewa-T surface in Freon-ll3. The obser ea 
behavior of this process is analyzed in this section by 


studying the different sections of the boiling curve. 
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Gewa-T Surface 
in Freon-113 


aeenetotemt Bolling Point 


boa h 


B: Natural Convection 

-C: Mixed Boiling (Activation of 
Nucleation Sites) 

C-D: Active Boiling 

D-E: High Flux Nucleate Boiling 
E-F: Stable Nucleate Boiling 


Piguse™ 7. Typical Nucleate Pool-Boiling Curve 
for a Gewa-T Surface 
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From point A to point B, a continuous increase Gnegete 
flux, accompanied by an increase in wall superheat (T. — Toad? 
is observed when heat flux iS increased. This region corres- 
ponds to natural-convection process, as also verified dumm 
experimental runs with the absence of bubbles along both 
enhanced surface and unenhanced ends. 

From point B to point C, a drastic reduction wine. 
superheat was observed while the heat flux continuously 
increased. This region was somewhat unstable as more and 
more nucleation sites became active. This region is known 
as the mixed boiling region, where transition from natural- 
convective boiling to nucleate pool boiling takes place. At 
these heat fluxes, visual observations revealed that the 
bubble formation was quite random. While some portions of 
the tube showed bubbles, other portions showed no bubble 
formation. The wall superheat continued to decrease until 
all nucleation sites became active, as indicated by point D. 

Once All available nucleation sites have been activated, 
the wall superheat starts to increase with increasing heat 
flux as shown by region D to E. When boiling in this region, 
the entire tube showed very high density of bubble formation. 
At this point, it is important to point out that aGuringe eae 
experimental runs, for the entire range of power settings (7.5to 
485 W), little or no nucleate pool boiling on the unenhanced 
surface of the ends was observed (see Figures 8 and 9) for 


the different test sections tested. This oebservation 
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Photograph of Gewa-T Surface at High Heat Flux 


8. 


Figure 
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confirmed the validity of the mathematical model assumed for 
Geta ceducticon. However, a small rate of bubble formation 
was evident at the tip of the unenhanced fin end; these 
bubbles were in fact generated at the base of the fin (i1.e., 
between the Teflon plug and the inner surface of the fin). 

The vapor bubbles generated in the inter-fin cavities 
gage Carried up into the channels around the tube due to 
buoyancy forces. As these bubbles travel up, they activate 
more and more nucleation sites as well as they help for the 
removal of larger bubbles from the channels. This phenomenon 
is generally referred to as the "Chimney Effect." 

If the heat flux is too large (see Figure 7), the rate of 
vapor bubble generation may be too large that the interfin 
channels could be full of a vapor blanket that may disrupt 
the flow of liquid to the hot metal surface where active 
boiling takes places. At this point, further increase in 
heat flux may result in excessive wall temperatures. Under 
these conditions, vapor bubbles were seen ejecting radially 
outward from the test surface. Line D-E represents the boil- 
ing regime just described. It is clear that the slope (i.e., 
heat-transfer coefficient) at E is considerably smaller than 
that at point D. This decrease in heat-transfer coefficient 
is assumed to be due to the existence of an insulating vapor 
blanket. 

When heat flux is gradually decreased, the curve follows 
a different path (D-F) as shown in Figure 7. This is due 
to the existence of stable nucleation sites that will remain 


active for a wide range of heat fluxes. 
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E. REPEATABILITY WHT VPRey Tous ss 

Following modifications made, especially to the data 
gathering means (1.e., automatic data-acquiSition systemeuame 
newly-written computer programs), it was necessary to take 
data for tube 1 (this is the same tube tested by Hernandez) 
to check the present data with Hernandez, data under similar 
conditions. Even “though vameec-e. m= different runs were 
made during this stage through debugging of programs and 
system familiarization, only run number 5 will be presented 
in this thesis with a comparison to Hernandez data. 

Figure 10 shows the comparison between the present data 
and the data of Hernandez. A very good agreement is evident 
as the two sets disagree by less than 2 percent for the 
nucleate boiling regime. However, considerable disagreement 
exists in the region of increasing heat flux before total 
nucleate boiling. This is believed to be due to the rather 
random behavior of the boiling process during this regime 
aS it will be presented in the next section, even two consecu- 
tive runs for a given tube would not give the same results 


during the unstable boiling regime. 


F. REPEATABILITY STUDY OF PRESENT Serr 

Figure ll shows data from two different runs completed 
on two different days. As discussed earlier, a considerable 
disagreement can be seen before the tube undergoes complete 
nucleate-boiling process. However, for the complete 


nucleate-boiling process, the data from different data runs 
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Pig mcomwititmyO.2 percent at a heat flux of 30,000 to 

moO, 000 W/m? . These data runs show considerably good repeata- 
memiyeOf Enis experiment. (Similar agreements were also 
feund £Or all other tubes.) 

Also, these data show a considerable enhancement of 
Gewa-T tube performance over the "Smooth" tube data of 
Hernandez. For example, at AT = 3.1 K, the Gewa-T tube 
gives a 4.65 times larger heat flux than the smooth tube. 

As also discussed by Lepere [Ref. 9] and Hernandez [Ref. 
8] this enhancement in boiling performance is mainly attri- 
buted to the existence of re-entrant cavities in the Gewa-T 


tube. 


G. BOILING PERFORMANCE OF SMOOTH TUBE 

To compare the performance of Gewa-T tubes to smooth-tube 
performance, data of Hernandez were chosen. It 1S important 
ae point out that his smooth tube did not consist of a commer- 
cially-available surface. Instead, Hernandez obtained his 
"Smooth" surface by machining fins off of a Gewa-T tube. 
He did not remove fins completely; the remaining portions on 
mae tube measured about 0.1 mm. 

Figure 12 shows Hernandez'sdata (heat flux as a function 
of wall superheat) for only decreasing heat fluxes. A 
least-squares-fit curve generated for these data points will 
be used in future figures representing data on tube 1. Since 
tubes 2 through 7 have a larger diameter, it was necessary 


memrevise Hernandez's data to include the difference in 
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Iareter. "Hee this puroese, a Semi-empirical correlation 


Bevelooued by Cornwell et al. [Ref. 11] 1s used: 
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Using this relationship, Hernandez's data were revised 
and a least-squares-fit curve was generated. This curve is 


gees plotted in Figure l2. 


peeeOoOLLING PERFORMANCE OF GEWA-T TUBE WITH SHROUDS 
iee Same four shrouds tested by Hernandez were used 
during this investigation. These shrouds had the following 


Pecmecry : 


a3 


Shroud Number Upper Opening Lower Opening 


al Sie 60° 
Z 60° Sun 
3 30s BO 
+ 60 Soe 


All of these shrouds had been designed for tube 1; thus, 
no shrouds were used for tubes 2 through 7. 

Figures 13, 14, 15, and 16 show data on tube 1 with Shea 
numbers 1, 2, 3 and 4, respectively. As also reported by 
Hernandez {Ref. 8], all the shrouds improved boiling per- 
formance at low heat fluxes compared to plain Gewa-T tube. 
Table 2 lists performance calculations for tube 1 with various 
shrouds. It can be seen that the 60° x8.5° shroud (number 4) 
produced the maximum improvement of 253 percent in the boiling 
coefficient over smooth tube at a heat flux of 10,000 W/m 
At this heat flux, the 30° x30° shroud resulted inva 
percent improvement in the boiling heat-transfer coefficient. 
At a heat flux of 50,000 mee all shrouds resulted in about 
190 percent improvement, except the 30° x 30° shroud which 
resulted in a 180 percent improvement. 

At very high heat fluxes, for example at q = 100,000 W/mee 
the 60° x60° shroud showed a maximum performance increase of 
106 percent; the 60° x 30° shroud resulted in a 49 pepcere 
improvement; the 60° x8.5° shroud resulted in a le eerece. 


increase; and the 30° x 30° shroud resulted in an 13 percent 
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improvement. It can be seen that the values just presented 
are somewhat proportional to the summation of the two shroud 
angles. This behavior indicates that at very high heat 
fluxes, too small a shroud opening can obstruct the outflow 
Of vapor and/or Ehewiaela wets tiapace 

However, it must be noted that at all practical heat 
fiuxes (10,000 to cunmous ae all the shrouds resulted in 
a considerable enhancement in the boiling performance. The 
60° x8.5° shroud appears to be the best selection for all 


practical heat fluxes. 


Jig PERFORMANCE OF GEWA-T TUBE WITH WIRE WRAPS IN INTER-FIN GAP 
A series of runs were made with two to five 0.1 mm-diameter, 
copper wire wraps in the inter-fin gap of tube 1. Figures 
17, 18, 19 and 20 show data with wire wraps of 2, 3, 4 and 
5, respectively. It can be seen that all of these curves 
show increases in boiling performance at all wall superheats 
for decreasing heat flux values. This enhancement is most 
probably due to the increase in the nucleation sites within 
the inter-fin cavities. However, there is a small reduction 
in performance in the natural convection process probably 
owing to that wires could act as an insulating material due 
to the thermal contact resistance between the wires and the 
inter-fin area. 
Table 2 also shows that the use of three wires resulted 
in the best performance increase. This combination resulted 


in a 341 percent improvement in the boiling heat-transfer 
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coefficient over the smooth tube value at a wall superheat 
oe 1.4 K. 
J. BOILING PERFORMANCE OF TUBE 1 WITH 60° x 8.5° SHROUD AND 

THREE WIRES 

Data were also taken to combine the best shroud and the 
Bese Number of wire wraps in the inter-fin channel. For 
this purpose, three wire wraps were added to the inter-fin 
channel and data were taken with the 60° x8.5° shroud. As 
Figure 21 and Table 2 show, this combination resulted ina 
performance yet superior to the use of three wires or the use 
of 60° x 8.5° shroud separately. Table 2 shows that this 
arrangement resulted in 481 percent increase in the boiling 
heat-transfer coefficient over the smooth-tube value at a 
heat flux of 10,000 ce 
eer PECTS OF FIN-TIP GAP AND FIN DENSITY OF GEWA-T 

/ TUBES ON BOILING PERFORMANCE 

As described earlier, six Gewa-T surfaces were tested 
during this stage of the investigation; tubes 2, 3, and 4 
imeeea fin density of 0.75 fins/mm, and tubes 5, 6 and 7 
with a fin density of 1.02 fins/mm, and each set of tubes 
feeeierin—tip gaps of 0.15, 0.25 and 0.33 (see Table l). 

Micrographs were obtained for all enhanced surfaces on 
meeeeanning Electron Microscope with a magnification of 50X. 
Figure 5 shows micrographs representing tubes 2 and 5. 
Computed surface areas and area ratios (with respect to 


smooth tube) are listed in Table 3. 
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TABLE 3 


Computed Surface Areas and Area Ratios 


Tube Cavity and Fin. Total Surface Area AD/AsmoorTH 
Surface Area (m ) (me) 
1 pag aT x10” mene Lon eae 
2 288.01 x10° 10.944 x10 2.74 
3 OEE lol - mieose 0 ~ 2.91 
4 285.07 x10 ° MOMceenx 10° - Day 
: 256.04 x10 ° Meena. 10) 3.60 
6 GION 12.998 Salo Peo 
7 eel x 10m ems x 10m. oe 


It can be seen that tubes with a higher fin density 
(1.02 fins/mm) result in larger area ratios than the tubes 
with lower fin density (0.75 fins/mm). However, for a given 
Bimecaensity, the variation of the area ratio with fin-tip 
gap is not conclusive. 

Figures 22 through 27 show data taken on tubes 2 through 
7, respectively. These results are also listed in Table 2. 
This table shows that tube 6 (with a fin density of 1.02 
fins/imm and a fin-tip gap of 0.25 mm) resulted in the highest 
heat flux at all values of wall superheat. 

Momim@vestigate the effects of fin density and fin-tip 
gap on boiling heat-transfer performance, the results pre- 


Pemtea 1n Figures 22 through 27 (which showed variations of 
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heat flux as a function of wall superheat) were re-plerr-e 

on a different basis as shown in Figure 28. This figure shows 
the variations of the enhancement of heat-transfer coefficient, 
over the smooth-tube value, with wall superheat for all six 
tubes. It is evident that tube 6, which has a £in densi, 

of 1.02 fins/mm and a fin-tip gap of 0.25 mm results inveee 
maximum enhancement at all wall superheat values. 

Tube 5 shows poorer performance than tube 6. This poorer 
performance is believed to be due to the fact that the smaller 
fin-tip gap disrupts the flow of liquid into the channel 
and/for the vapor flow out of the channels. Tube 7 (which has 
the largest fin-tip gap of 0.35 mm) shows the poorest pea. 
formance among the set of these tubes with a fin density Ge 
1.02 fins/mm. As visually observed during experimental 
runs, this tube ejected a considerable portion of vapor 
bubbles radially along tube periphery. This = orem. 
departure of vapor bubbles results in a lower performance 
by producing a thicker liquid film at the hoe surracee 

In comparison with the results for tubes 5, 6 and 7, the 
curves representing tubes 2, 3 and 4 show that for a given 
fFin-tip gap, the performance has decreased with decrease in 
fin density. This observation is believed to be due to the 
presence of a greater number of channels (for a given length 
of tube) in the higher-fin-density tubes. Again, for a fin 
density of 0.75 fins/mm, the 0.25 mm £in-tip gap (eubon 


resulted in the best boiling performance. The same 
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explanations given above should be valid for the poorer per- 
formance of tubes 2 and 4. 

It can be seen that all tubes show smaller enhancement 
ratios at low wall superheats, and the above-mentioned 
differences in performance are more pronounced. On the 
other hand, all the tubes show greater enhancements at larger 
wall superheats and they seem to approach a common enhance- 
ment ratio. The increase in enhancement is believed to be 
due to the increased chimney effect, which helps to remove 
vapor bubbles from channels resulting in a thinner liquid 
film at the hot boiling surface. When heat flux is further 
increased, vapor can provide a blanket around the hot 
surface, which may eventually result in a critical heat 
flux. It appears that the critical heat flux iS not stroemame 
dependent on the fin density or fin-tip gap. 

The results shown in Figure 28 were also replotted as 
shown in Figure 29 to further study the performance of these 
Six tubes. This figure shows the variation in enhancement 
ratio obtained beyond the area increase (1.e., the ratio of 
heat-transfer enhancement/area ratio) as a function of wall 
Superheat. It is evident that obviously all tubes resulted 
in enhancement ratios beyond the area increase. 

In summary, a fin density of 1.02 fins/mm resulted ig 
performance superior to a fin density of 0.75 fins/mm. Also, 
in both cases, a fin-tip gap of 0.25 mm resulted inetieeeeee 


performance. 
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Vo CONCLUSJGNS 


(l) Experimental observations revealed a non-uniform 
temperature distribution along test section. The bubble 
generation was seen to start at the two ends and migrate into 
the center as the heat flux was increased. No clear expienee 
tion is possible for this temperature distribution. 

(2) Results of the present investigation showed very 
good repeatability with the results of Hernandez obtained 
under similar conditions. 

(3) The use of shrouds improved the boiling heat-transfer 
coefficient of Gewa-T surface at all practical heat fluxes? 
However, a considerable reduction in heat-transfer coefficient 
was seen at very high heat fluxes (> 100,000 w/o Among 
the four shrouds tested, the shroud with 60° x8.5° openings 
resulted in the best performance verifying Hernandez's 
results. This shroud resulted in a 253 percent increase in 
boiling heat-transfer coefficient, over smooth tube, at a 
wall superheat of 1.8 K. 

(4) The use of 0.1 mm-diameter wires (2 to 5) improved 
boiling heat-transfer coefficient at all practical heat 
fluxes. The use of 3 wires resulted in the best performance 
with a 341 percent increase in the heat-transfer coefficient 
(over the smooth-tube value) at a wall sSuperheat of 1.4 K. 

(5) Further improvement in boiling performance was 


obtained by using the 60° x8.5° shroud with three wire wraps. 
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This combination resulted ina 481 percent increase in boil- 
ing heat-transfer coefficient at a wall superheat of 1.1 K. 

(6) The use of wires appears to be superior to the use 
of more expensive shrouds. 

ee Menaqmene fan—tap gaps of 0.15, 0.25 and 0.35 mm 
tested, the 0.25 mm gap resulted in the best boiling 
fear Oormance . 

(8) The Gewa-T tubes with a fin density of 1.02 fins/mm 
performed superior to a fin density of 0.75 fins/mm. Ata 
wall superheat of 3.5 K, the lower-fin-density tube increased 
the boiling heat-transfer coefficient by 80 percent, while 


the higher-fin-density tube increased by 103 percent. 


Te, 


VI. RECOMMENDATIONS 


(l) Investigate the reason for the existence of a 
longitudinal temperature profile and take) appeepeiae- 
steps to alleviate this problem. 

(2) Study the possibility of adding roughness aayeie 
inter-fin cavity of Gewa-T tubes.” This operation sneum 


be performed before the rolling process to obtain the T- 


shaped fins. 


(3) Repeat runs made during this investigation with 


water as the boiling medium. 
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APPENDIX A 


DATA REDUCTION PROGRAM 


eee ene ORS 


DATE: May °3. 1984 
REVISED: Suiy 20. !9a4 
REEP 


PRIMER 1S 7 
PRINT USING “4X.""Select option: ?'" 


PRINT USING “6X."°"0 Taking data or re-processing previous data” 


PRE USING (BX. I Plotting data ~” 

INPUT Idp 

IF Idp=0 THEN CALL Main 

IF Idp=1 THEN CALL Plot 

END 

SUB Main 

ColeGe7 wl C7) 

DIM Emf(9),.TC(9) 

DATA 0.10086091 .25727 .943569, -757345.38295, 78025595.31 
DATA -9247486589.6.97688E+11,-2.66192E+13.3.94078E+14 
READ C(«) 

BEEP 

(evi soeeLec) (UBE UD <O0=0:.7Ssin, 1=)20 ind”. ltd 
peeea=j0) THEN 
D1 =.0159 
D2=.0133 
Diao Ss 
Do=.0190 


Diameter at thermocouple positions 

Diameter oreatest section to the base of fins 
Inside diameter of ‘tnenhanced ends 

Nutside diameter of unenhanced ends 

Length of enhanced surface 

Kength of unenhanced surface at the ends 


ne 9 ome me 


D1=.01786 

P2=20231 

Di=7 02032 

Do=.02222 

L=,0508 

Lt 303175 

END IF 

Ee-=hur(0o-01)/2 

Asoc pe 2-01 2)/4 

P=PT<Do 

Keu=385 ' Thermal conductivity of Copper 
PRINTER IS 701 

Clana, 09 

BEEP 

ive ENTER MONTH. DATE GND TIME «MM: DD:HH:MN:SS)” .DatesS 
Dewees 10 sDates 

BUIPOie7Z tas 1D” 

ENTER 709:Date$ 

eon) ee Honcunwe alate Ancdmrimnem = «9 Dates 
PRINT 

eee BSrnGe 0k. “NOTE: Program name: BUIL"’” 


INPUT “ENTER DISK NUMBER”. Dn 


cae USING “16x, “Disk number = "".ZZ7"°:Dn 

INPUT “ENTER INPUT MODE (N=2054A,!=FILE)".Im 
Im=Im+1 

TF Tm=!t THEN 

BEEP 

INPUT "GIVE A NAME FOR THE RAW DATA FILE’.D *1les 
ProtieUoeiGs 16x, “New ile name: “".14A":D_ files 
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SEEL 
INPUT “ENTER FLUID CO=FREQN=1 1301 -NAle ee 
CheCi + 


\ 
GUTPUT wFile:Dates.c€ 
Boe 

BEEP 

yea “GIVE THE NAME OF THE EXISTING DATA FILE” .D_files 
BEEP 

The=8} 

INPUT “ENTER [t IF USING FILE NUMBER < 25°. [he 

ey USING “16X.°"Old file name: ”° 214A 340m eles 


INPUT “ENTER THE NUMBER OF RUNS STORED’ .Nrun 

ASSIGN @File TO D_fileds 

ENTER #File:DoldsS,Cf 

PRINT USING “16X%.°"°This output for ftle created at "FI4A Guo 
END IF 

BEEP 

INPUT “GIVE A NAME FOR PLOT FILE’.P_fileSs 

CREATE BDAT P_files.5 

ASSIGN *Plot TO P_files 


IF Cf=1 THEN PRINT USING "16X%."°"F luid : Freen= liga 
IF €f22 THEN PRIN USING 16xX.. Pie : Watecue= 
BEBE 


INPUT “ENTER BUTPUT VERSION ¢1=SHQRT.0=LUNG)”.Iov 
IF Iov=0 THEN Iov=2 


BEEP 

INPUT “ENTER NUMBER OF DEFFEGIIVE TCS “0=DEeAULD Sia 
IF Idtc>O THEN 

BEEP 

[Poicte= es iretl ; 

INPUT “ENTER DESFEGI IVE FC Buca sear 

Oe "“16X." "TC rs deffective at location ~~" 7ueeecmel 
NM 

IF Idtc=2 THEN 

INPUT “ENTER DEFFECTIVE TC LOCATIONS” .Ldtc!.Ldtc2 

PRINT USING "16X.""TC are deffective at locations “”.D,.4X.D° :Sate team 
END IF 

END IF 

IF Im=! THEN 

BEBE 

INPUT "ENTER BAROMETRIC PRESSURE (1nHg)"’.Patm 
Patm=Patm+2=1500/13500 

QDUTPUT OF ite;Patm 

Bee 

ENTER 3File:Patm 

IF Ihc=1 THEN Patm=Patm+2*1500/13500 

END IF 

Patm=Patm*25.4*14.696/760 

Tsat=FNisat¢(Cf ,Patm?> 

PRINT USING “!16X.°"Saturation temperature = “",3D.DD.°" (Deq G)ueems 
j= 

Sx=0 

Sy =0 

Sxs=) 

Sxy =) 

pee 

tNPUT “HIT CONTINUE WHEN READY”. Ok 
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e200 Repeat. — 


Ee Ab 
eco 
E30 


2240 
225) 
2260 
eo 
2280 
Ze) 
2300 
20 
2320 
eo) 
2340 
23350 
2360 
27/0 
2380 
2390 
2400 
2410) 
, 2420 
2430 
2440 
2450 


a0 
252) 
aU 
2940 
2500 
2560 
B00 
Zo00 
230 
2600 
2610 
2620 
2630 
2640 
2650 
2660 
26/70 
2680 
2690 
2700 
a 10 
e720 
eto 
2750 
eT U 
2760 
eT 0 
2780 
2730 
2800 


Ceo Ga) THES! 

DUTPUT 709:"98R AF29 ALSB VRS” 

FOR [=6 TO 3 

Hour e7og- "AS SA” 

Suun=0 

FOR Ji=1 TO 20 

FNTER 709:E 

Sum=Sumte 

NEXT Ji 

Emf (I) =Sum/20 

NEXT I 

Bea 

INPUT “ENTER HEATER VOLTAGE. Vh (Y)" Vb 
BEEP 

INPUT "ENTER RESISTOR VOLTAGE. Vr (V)" Ur 


oO 
eee oes pnt (=, 75. 7h 
END IF 
Twa=4 
FOR I=0 10 9 
IF Idtc>0O THEN 
het =_Ldtc!l OR «<141)=Ldte2-THEN 
el ==9.99 
GOTO 2490 
2/0) 103 
END IF 
TCI) =FNTIvsviEmt (I)? 
Pael<S IHREN S lwa=lwatt Cl) 
NEXT I 
Twa=Twa/(8-Idtc) 
T1ld=T(9) 
Tv=1(8) 
Rho=FNRho(Cf .Tsat) 
Mu=FNMu(Cf,. Tsat) 
K=FNK (CE, Tsat) 
Ga-=- NCoCC#H , Isat) 
Beta=FNBeta(Cf.Tsat) 
Ni=Mu/Rho 
Alpha=K/(Rho«Cp) 
Pr=Ni/Alpha 
Amp=Yr/2.031 
Q=Vh*Amp 
Tw=Twa-U*LOG(D2/D1)/(02*PI*KcueL) 
metab = tu isat 
IF Thetapb<0 THEN 
bee 
DISP "TWALL<TSAT (CHIT CONTINUE WHEN READY)" 
PAUSE 
GOTO 2210 
END IF 
Hbar= 190 
Pe=(hbar=P/(KcusA)) .5Le 
Tanh=FNTanh(Fe) 
Tneta=lhetab«Tann/Fe 
Xx=(9.831*BetaeThetab«Do 3*Tanh/(FeeNi¢Alpha)) .156667 
Yy=(1+¢€.559/Pr) (9/16)) (8/27) 
Hbarc=K/Do*( .6+.38/7<=Xx/Yy) 2 
IF ABS( (Hbar-Hbarc)/Hbarc)>.001 THEN 
Hbar=(Hbar+Hbarc)<.5 
GOTO 2720 
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2 0) feviDe 2 

2320) Ql=(Hbar-P<K cur) .5-Therab<Tann 

2630 Qc=0-2</)) 

2340 As=rPl=-D2-L 

2850 Qap=Qc/As 

2360 IF Iov=i THEN PRINT USING "10X,0D.2X,3D.DD.2X%,.2(Z.4DE.2x) DD DD “eee 
PT Udpmrne tao 

3/70 = 6oTF 6Lovw=2 THEN 


2380 PRINT USING “70K.” Data set mumper -= 7) 0ee 
28390 pod USING “10K. Poei tien  numoer a 1 ze S 4 5 b 
7 eeeoes 
2900 i) vous "“10X.°°’ Temperature : *" 8(DD.DD.1X)": 700). 1 10 G2 
CAD C57 wey) 
2310 ata ater “TOX, oT lita Tvapr Twa rite! Qdp The tain aaa 


2920 PRINT USING “SX 3°30. DD 2x0 (202 405 2 eevee ila) Tv. Twa. Hbar ,.Qdp. Inhetab 
23300 “END IF 

2940 IF Im=1 THEN 

E20 wise ee 

2960 INPUT “OK TO STORE THIS DATA SET ¢1=Y.0=N)?” Ok 
2970 END IF 

23980 IF Ok=t OR Im=2 THEN J=J+! 

2990 IF Ok=1 AND Im=1 THEN OUTPUT ®&File:Emt (C=) Ver oVa 
3000 IF Im=2 OR Ok=1 TREN DUIPUT SPlot:@dp. fhetab 

3010 IF Im=1 THEN 

Su ZU ye peee 

3030 INPUT “WILL THERE BE ANOTHER RUN (1=Y,N=N)?",Go_9n 
3040 Nrun=J 

3US0 IF Go_on<>! THEN 3119 

3060 WAIT 300 

3070 IF Go_on=! THEN Repeat 


3080 ELSE 

3090 IF J<Nrun+i THEN Repeat 

3100 END I] 

3110 IF Im=! THEN 

3120 BEEP 

3130 PRINT F 

3140 PRINT USING “!0X.°"NOTE: °".ZZ.°° data cuns were stored in fale see 
120 Files 

2750) VEND IF 

3160 BEEP 

3170 PRINT 

3180 PRINT USING “10X.°"NOTE: “".2Z."" X-Y vairs were stored in Blo usa 
fo Oh I aces le> 

3190 ASSIGN 4File TO < 

2200 ASSIGN @Plot TO + 


Be) eee 
3220 <INPUT “LIKE TQ PLOT DATA ¢1=7,0=N)?", 0k 
Ae IF Dk=1 THEN 
240 CALL Plot 
3250 END IF 
3360 SUBEND 
3270 DEF FNMutCf.T) 
3280 IF Cf=!1 THEN 
53230! RANGE: 26 TES 0sDecee 
3300 Mu=8,.9629819E-4-T<¢€1,.1994609E£-S-T+5.566829E-8) 
3370 END IF 
3320 IF Cf=2 THEN 
3330 A= 247, SiC is) 3a, 5D 
3240 Mu=2.4E-5<!10 4 
5350 END. IF 
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RETURN Mu 

FNEND 

Dep oNCe tls, |) 

IF Cf=1 THEN 

ReENGe. 0 60 360 DEG C 


es) 3-1 + 1-05 eee -4+1 . 7207792E-6-T) 


END IF 

oe TEN 
Pere ec ee oem om cob FES he. 7 1S20E Sy 
EN F 

RETURN Cp#1000 


FNEND 
DEF HENRne< Cf . 7) 
fe eC le HEN 
PeawNGe sou 10 80 DEG C 


Booed a7 Ge too 2, clopodotii=2.35/6cJ71E-3) 
END IF 

[IF Cf=2 THEN 

Ro=999 .52946+T*¢( .01269-T<¢(5. 4825136 -3-T<«1.234147E-5)) 
END IF 

RETURN Ro 

FNEND 

CEP ENP TC, ID 

Pr=FNCo(Cf .T)=FNMuC(Cf.T>/ENKCCE.T) 

RETURN Pr 

FNEND 

DEF FNK(CF,T) 

IF Cf=1 THEN 


ANGE 30 10° 30 DEG C 


epee idsesee-2 -1= 2.22 14206E-4e1~2.. 38095245 -3) 
ele 

IF Cf*2 THEN 

Kalter. boovers.15 

ee tae Oo IO ( 1. S007 -X*¢ .52577-.0/7344"X))) 


- END IF 


RETURN K 

FNEND 

DEF FNTanh¢X) 
P=EXP(X) 

Q=EXP(-X) 
Tanh=(P-Q)/(P+Q) 
RETURN Tanh 

FNEND 

DEF FNTvsv(V) 

CO wte/  'CC7) 
T=C(0) 

FoR a ho? 
T=T+eC] =) I 

NEXT I 
T=71+8.94643E-1-T<(2,4547896E-2-T«1.3722583E-4) 
RETURN 7 

FNEND 

DEF FNBeta(Cf.7) 
Rop=FNRho(Cf ,T+.1) 
Rom=FNRho (Cf ,T-.1) 
Beta=-2/(RoptRom)*(Rop-Rom)/.2 
RETURN Beta 

FNEND 
VEPOENIsatcCt VP) 
ft) 

Th=110 
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So7 00 s= Ol ehee 
3980 Pe=FNPsar(Cf.Ta) 

3990 IF ABS CCP-Pc)/P)>.0001 THEN 

4000 IF Pe>P THEN Th=Ta 

4010 IF Pc<P THEN Tl=Ta 

4026 G18 3370 

ANS05 END ie 

4040 RETURN Ta 

4050 FNEND 

4060 DEF FNPsat(Cf.Teo) 

406%) If €f=1- THEN 

4070 T=To~!.8+32+459.6 

4080 P=10°(33.0655-4330.38/7T-9.2635*LGT(T)+2.0539E-3*T) 
4081 END IF 

4082 IF Cf=2 THEN 

4083 P=Pvstuw(Tc) 

4085 END IF 

4090 RETURN P 

4100 FNEND 

4110 SUB Plot 

A120 De Gap 

4130 INTEGER [1 

4140 PRINTER IS 1 

4iS0- -tdvet 

4160 BEEP 

4170 INPUT “LIKE DEFAULT VALUES FOR PLOT (1=#Y.02N)?",Idv 
4180 BEEP 

4At30 PRINT USING “4x. “Selisetlptrten: ) 

4200 -PRINT USING “4X9 “O "aeversicodelta-i 

4210 ‘PRINT USING “4X. ee sieverscus delra-i. 

4220 PRINT USING “Gke* 2." versus a” 

4230 PRINT USING “4XK.°"3 h=ratio versus delta-i 
4240 INPUT Opo 

4250 IF Opo=3 THEN 


4250 Beer 
4270 INPUT “SELECT TUBE DIAMETER (0= 7 Sie) sierra 
4280 END IF 


4290 PRINTER IS 705 

4300 IF Idv<>1 THEN 

430) seer 

4320 INPUT “ENTER NUMBER OF CYCLES FOR X-AXIS™.Cx 
4330 BEEP 

4340 INPUT “ENTER NUHBER GF CYCLES FOR Y-ex ion, 
4350 BEEP 

4360 INPUT “ENTER MIN X-VALUE CNULTIPLE OF {0)9"°,Xmin 
So Ceee 

4280 INPUT “ENTER MIN Y-VALUE (MULTIPLE OF [9)".Ymin 
41690) SESE 

4400 Cy=3 

4410 IF Opo=0 THEN 

4420 Cxr=2 

4430 Xmin=l 

4449 Ymin=!000 

4450 END IF 

4460 IF YDpo=! THEN 

4a70 .Cx=2 

4480 Xmin=} 

4490 Ymin=100 

4500) ENO IF 

4510 IF Dpo=2 THEN 

AD 20) n= 3 
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4530 Xmin='000 

4540 Ymin=!00 

A550 END IF 

4560 IF Ypo=3 THEN 

Hoy 0) sve x= 2 

4580 Cy=3 

4590 Xmin=! 

4600 Ymin=1 

4610 END IF 

4620 END IF 

4630 BEEP 

4640 PRINT “IN:SP1;:IP 2300.1900,3300,6800;" 
sole iNt ‘SC 0.100,0.100;TL 2,0:" 

4660 Sfx=100/Cx 

4670 Sfy=100/Cy 

4680 PRINT “PU 0,0 PD" 

4690 Nn=9 

ZOO FUR f=) 18 CxFl 

wat Kat=kmin*x10° (1-1) 

4720 IF I=Cx+1 THEN Nn=! 

4730 FOR J=1 TO Nn 

4740 IF J=1 THEN PRINT "TL 2 0” 

A750 8ir Jee THEN PRINT “TL 71 0” 

4760 Xa=Xat+J 

S720) X=GT(Xa/Xmin)*St x 
4780 PRINT "“PA':X."°.0: XT:" 
4790 NEXT J 

4300 NEXT I 

4310 PRINT “PA 100.0:PU:” 
4820 PRINT “PU PA 0,0 PD" 
4830 Nn=9 

“oa; FOR i=! 1G Cy+t 

Heo Yat=1min=\0 €i-1>) 
4860 IF I=Cy+1! THEN Nn=! 
4870 FOR J=1 70 Nn 

4380 IF J=! THEN PRINT “TL 
eee0 if J=2 THEN PRINT “IL 
AV Ya=lat=J 

4910 Y=LGT(Ya/Ymin)*Sfy 
BS20 8 PRINT “PA O0.°:Y, YT” 
4930 NEXT J 

43540 NEXT I 

4950 PRINT "PA 0,100 TL 0 2" 

4960 Nn=9 

4970 FOR I=! TO Cx+l 

foo) Kat-amim=10 ({[-1> 

e250) (F f=Cx-! THEN Nn=] 

S000) *-FOR J=1 10 Nn 

ete =) THES PRINT “TL O 2” 

Bec IF J>1 THEN PRINT “TL O 1" 

5930 Xa=Xat*J 

5040 xX=LGT(Xa/Xmin) «Sfx 

Sot rRINT “PAX. °. 100: XT” 

5060 NEXT J 

S070 NEXT I 

SuSUmeeRIN, “Pa 100,100 PU PA 100.0 FD" 
5090 Nn=9 

5'00 FOR I=! TO Cy+! 

Si Yat=Ymine«10° <«(1-1> 

5120 IF T=Cy+! THEN Nn=! 

5130 FOR J=1 TO Nn 


2 
1 


0 
QO" 
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iF J=l 7HEN PR Flee. 
IF J>1 THEN - PRINTS Teer 
Ya=YatrJ 
Y=eLGTCYa/Ymin)«Sfy 

PRINT “2D PA SOON oye) i 
NEXT J 

NEXT f 

PRINT “PA 100.100 PU" 
PRINT “PA Us] 2 Sh ies 
Li =EG I <min 

FOR I=1 TO Cx 
Xa=Xmin*«l]0° (I-11) 
X=LGT(Xa/Xmin)*Sf x 

PRINT “PAs ae 

IF T1>=0 THEN PRINT “CP —2,-2:CB10°PR —2.2 oe 
IF [i<0 THEN PRINT “CP —-2.-2°-LBI0:FR 0.2 26 eae 
Lrata! 

NEXT 1 

PRINT "PU PA 0,0" 
T1=LGT¢¥min) 

Y1o0=10 

FOR I=1 TO Cyt 
Ya=Ymine*10 > (1-1) 

Y= Re oa eee 

PRINT “PA -05 

PRINT “CP -4,- os. EB OTeR <2. 2 08 2 bia 
Pee | 

NEXT I 

IF Idv<>!t THEN 

BEEE 

INPUT ENTER X-Eas6e = Glabe > 
BEEP 

INPUT “ENTER Y-LABEL”,.YlapbelS 
ELSE 

IF Opo=0 THEN 
XlabelS=" Tw lsat) Ck) 
YlabelS=""q (W/m 2)" 

END IF 

TF ODpo=! THEN 
XlabelS="Tw-Tsat ‘«(kK)?” 
YlabelS="h (H/m 2.8? © 

END IF 

IF Opa=2 THEN 

XlabelS="a CW7ym 2)” 
YlapelS="h -GH7 mick). 

END IF 

IF Bpo=3 THEN 
XlabelS="Tw-Tsat (K)” 
YlabelS="ht enh) /h( smooth)” 
END IF 

CND IF 


PRINT “SR 1.5.2:PU PA 50.-16 CP":-LEN(Xlabel3)7 27 ing Rely 
PRINT “PA -14.50 CP 0.":-LEN(YlabelS)/2*5/6: DI 071 72B i) leper ee 


Smut CP 0 2001p i 


Repeat: ! 


BEEP 

INPUT “WANT TO PLOT DATA FROM A FILE ¢!=Y.0=N)2" Dk 
IF Dk=1 THEN 

pee 

INPUT “ENTER TRE NAME OF THE DATA FILE” DO _files 
ASSIGN File TQ D_files 


8 8 


ayipel!) 
5760 
3770 
5780 
Syeshl, 
2800 


BEEP 
BEEP 
INPUT “ENTER THE BEGINNING RUN NUMBER" .Md 


BEEP 

Beis ENTER THE NUMBER SUF X=7¥ PAIRS STORED” .Npairs 
Beer 

PRiNteR LS 


PRINT JSiNG “4X, ""Seleet a symbol ; 

PRINT USING “4X.°"i Star a ee Sai ae 
PRINTeUSING “4X.°"S Cirele 4 Saquare’’’” 
PROENGMUGSING “4%, S  )Rombus 


PRINT SING “4X, °"°G “Rigmt-side-up triangle”’”” 
PRINT USING "“4X,""7 Up-side-down triangle”’’” 
INPUT Sym 

PRINTER TS 70S 

PRINT “PU DI" 

IF Sym=1 THEN PRINT "“SM*"’ 

IF Sym=2 THEN PRINT “SM+” 

IF Sym=3 THEN PRINT “Sto” 

IF Md>1 THEN 

FOR I=1 TO (Md-1) 

ENTEkwehile:Ya,Xa 

NEXT = 1 

Eb. LE 

pope UeNoatre 

ENIER 88h ple: va. Xa 

IF Gpo=1 THEN Ya=Ya/Xa 

IF Opo=2 THEN 


END IF 

Ie Gee=2  1HEN Ya=va/eENHsmootint Xa. ltd) 

X=LGT (Xa/Xmin)d«Sfx 

Y=LGT(Ya/Ymind=Sfy 

[eeovym> oe HEN eRINT “SM” 

tae oyM< om ineN PRINT) oR 1.4.2.4" 

PACU tie Cimon ha) 

Te -Sym>3 THEN PRINT “SR 1.2.1.6" 

IF Sym=4 THEN PRINT ees oe oa Se .e 0: 
IF Sym=5 THEN PRINT "“UCS3,0. 33, eo ore fone i: 
IF Sym=6 THEN PRINT "UCO.5.3.99.3.-8.-6. oo an ‘ 
IF Sym=?7 THEN PRINT “UCO,-5S. Ae Soe ooo oes 
NEXT 1 | 

Beare 

ASSIGN File TO < 

GOTS 5690 

END IF 

PRINT “PU SM" 

BEER 

INPUT “WANT TG PLOT A POLYNOMIAL (1=Y.0=N)?".Go_on 
IF Go_on=1 THEN 

Baer 

INPUT “ENTER LOWER AND UPPER X-LIMITS".X1IL.Xlu 
FOR Xx=0 TO Cx STEP Cx/200 

Ka=Xmin«l0 Xx 

IF Xa<Xll GR Xa>X!lu THEN 6380 

Ya=FNPoly (Xa) 

feo l yay imine ~oi ¥ 

X=(GT(Xa/Xmin) «Sfx 

LF Y<0O THEN Y=0 
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s36bU ir YO!1UU- “REN “oT soc 

6370 “PRINT "RPA 2X27 ene 

6380 NEXT <x 

6390 END IF 

6400 PRINT "PU PA 0.0 SPO" 

6410 SUBEND 

6420 DEF FNHsmooth(X.I td) 

6430 Hs=FNPoly(X)/X 

6440 IF Itd=1 THEN Hs=Hs*.33347 

6450 RETURN Hs 

6460 FNEND 

6470 DEF FNPoly(X) 

6480 Poly=-4.4123718E+2-X=(6.812391 7E+2-X#3./7416863E+2) 
6490 RETURN Poly 

6500 FNEND 

6510 DEF FNPvst¢Tsteam) 

6515 DIM K¢s) 

6520 DATA -7.691234564.-26.08023696.-!68.1706546 .64.23285504.-118.9646225 
6525 DATA 4.16711732.20.9750676.1.E9.6 

65530 READ K(*) 

6535 [=< lsteamt2/3.15)7647.. 

6540 Sum=0 

56545 FOR N=0 TO 4 

56550 Sum=Sumt+tK«(N)<«C1-T) (N+!) 

6555 NEXT N 

6560 Br=Sum/€T«€14+K (5) «€C1-T) +K (6) 401 -T) 2) - CF -TIS CK C7 ¥ CI -T) 2 +K 8d) 
6565 Pr=EXPC(Br) 

56570 P=22120000«Pr 

6575 RETURN P 

6580 °NEND 
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APPENDIX B 


SAMPLE CALCULATIONS 


Data Set Mumber fifteen of file GT6H29 corresponding to 


tube number 6 was chosen for the sample calculation. 


Pao LEST-—SECTION DIMENSIONS 


Cn)’. — 2 Seo Ome min) 
Hi - 130 SaaGun) 
fins/mm = 102 
iiter-fin gap = 0.25 Sum) 
D. =e SO mira) 
Dy = 17.86 (mm) 
D — nn 
O 
Di. = 20.32 (mm) 
1 
L., = 31.75 (mm) 
L = 50.80 (mm) 


B. MEASURED PARAMETERS 


V = 70.0 (volts) 
Vers = foe volts) 
R = Zo 0 Sle e(ohms) 
Ko =e Coe Ml.) 

T) =e ao CC) 

T. eer C) 

T 5 ane) So Lil oC) 


oul 


i =o 6s. 2 C) 


4 
T. = 55240" ( 7e) 
Te. = (535) 
TS =. “Oop te) 
T. = basso 4) eS) 
Toar by = “Siar (2 e) 


PROPERTIES OF FREON-113 AT SATURATIONS TE PER ei 


u = Amoden. tome (kg/m.s) 
0 = 1510.0 (kg/m>) 

eS = €973 (c/ko Ww 

k = 0.0715 (W/m.K) 


HEAT-EhUx -CALCULATICGN 


The heat transfer rate from the cartridge heater: 


Oe. = UVES) and I Vepo/R (amps) 
Te= 428/20031. =)" 27363 eames) 
Q., = 170.0 «2.363 = 6574 emis 


The cross-sectional area ©f& machined sends: 


= Z 2 Z 

Ac = 1/4 .0(Do D.) (m_ ) 
2 e -6 Z 
Aa = 7/4.0(0.02222 - 0.020325), = 63-4252) (m_) 


Se 


The tube outside wall perimeter: 


P = wD, = 3.1416x0.02222 = 6.981x10°° (m) 
The corrected length of ends: 
ee = ee 2 a ole 
Hey = Oroe7s + 0.02222 = 0.02032)/2.0 = 32.710" (m) 
Other thermophysical properties: 
= = -4 _ =); Z 
Peeeewyo = 4.948 x10 /1510.0 = 3.277 x10 (m/s) 
=o Z 
a = k/oC, =O Lo7 Glow. x 97s) = 4.866 x10 (m/s) 
, je Ter i (1509.77 ~ 1510.26) 
. ae 1510 AG. 61g = 4 7a 4k 
=O OOro2 (1/k) 
Pay ee ep 108 /4esec x10 ° = 6.734 
Temperature at the base of the T-shaped fins: 
E | 7 In(D,/D,) 
B AVG H Zee Ko 
8 
ee et 78-0) = 523.69 (°C) 
n=l 
_ 2 MO ly os 6) z 


oe 


0 =) a = 53.342 - 47.51 — eee 


Average heat-transfer coefficient ae Ehewunerinage- = 








ends: 
= yee 
5 hp V7 
g8D56,Tanh((; x ) Lol 
Cac 5 
aE — yi f2 
h 0: 60 -OREIER ss 
h = = ; ; 
De [1 + (0.5597 pe 
= 2 
h = 187.49 (W/m .K) 
Heat-transfer rate through unenhanced ends: 
Qn = nh PKA. Op Tanh mL 
a = ly ae 
ae oe 1/2 _ (287.489 x GE Ant = 1.608 x 102 
Cac 020715 «63.43 Pex aug 
Op, = 3.29 (W) 
Qross = 20. = 62,579 eine 


LOSS 
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Q = 165.41 - 6.579 = 158.83 (W) 

q = Q/A, 

An = si 223. male 50.8 lies = 3.687 malo oa 
q = —a2e:8s = 43.978 x10° (W/m*) 


Bros 7 x LO - 


The following are the results obtained from the computer 


by running the data reduction program (DRP) (see Appendix A): 


WAGE 10 Gan) 


QO — 
0, = Mero) Oli) 
2 
h = 187.39 (W/m .K) 


The small difference in results obey the inaccuracy of 
hand calculations with the use of about three significant 


eraits only. 
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APPENDIX 


UNCERTAINTY ANALYSIS 


The same data set (set No. 15 of file GT6H29) that was 
used for the sample calculation was chosen for the uncertainty 
analysis; therefore dimensions of test section, meaguned and 
calculated parameters found in sample calculation were used 
in this analysis. All uncertainties are presented as a 


percentage of the calculated parameter. 


A. UNCERTAINTY IN SOURCE HEAT=TRKANSPER- RATE 





O48 5= aT OW) I = Verc/R (amp) 
6V = SVepc =" 03 05> (volts) 
where: 
6 = uncertainty in measurement and calculation 
6I = SVppc/R = 0.05/2.031 = (0202353 (ame,) 
001 a2) Benny 2 
Quy = [(—) (ae 
V ‘le 
6Q 72 
Sl OO a 0 302 Sa _ 
70m = leat 70 ) (5 363) ] = L206 peErCenue 


SG 


UNCERTAINTY IN SURFACE AREA 


A, = 7D, L 6D. = oF = 0.1 (mm) 
6A, 2 D1 
- O.1 O.1 _ 
oar oc. * ‘song) } ence 


UNCERTAINTY IN AT 





Ls = Ty = Toa oT = OToag =e s (°C } 
SAT oT 2 ST oar 2272 
AT OFF —) + (= ml 
oe AT 
Bee 
SAT _ 0.15,2 0.15,2 - 
some 6. +.‘ ¢.18) | Wa oo 


UNCERTAINTY IN HEAT FLUX 


Gg = a 
ae 
er? p,. *8B, 2/7 
4 O77 2Qp Q17 2p = 
In this case, Qn, = 165.41 (W) and Qn = 3.29 (W) 
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= 2,70 ercemre om Qu: Assuming the same 


Therefore, Qn = 
Proportion in the uneecritaime, soe Qn! 


6 On = 0.02 6On = 0.02 x1.753” = WO? 
172Z 
car 7 Soe Bee < U2 02a _ 2 
q = l (seg BF) + ( 153-63? + (-0.0048) "] 
2d = 1.26 percent 


UNCERTAINTY IN BOILING HEAT-TRANSFPERBGOE? Plc las 


he 
bee = CSE 
: AEE 
SR (82 4 (- S8B)) %* 
"I q AT 
= eZ 
oa = ome eG) = + (20 osen = 3.65 percent 
h 
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nO. 


alle, 


2 . 
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